Scaling Theory for Migration-Driven Aggregate Growth 
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We give a comprehensive description for the irreversible growth of aggregates by migration from 
small to large aggregates. For a homogeneous rate K{i;j) at which monomers migrate from aggre- 
gates of size i to those of size j, that is, K{ai\ aj) ~ a^K{i\j), the mean aggregate size grows with 
time as i^/'^-A) A < 2. The aggregate size distribution exhibits distinct regimes of behavior 
that are controlled by the scaling properties of the migration rate from the smallest to the largest 
aggregates. Our theory applies to diverse phenomena such as the distribution of city populations, 
late stage coarsening of non-symmetric binary systems, and models for wealth exchange. 



Much attention has been devoted to understanding the 
irreversible growth of aggregates through binary coales- 
cence. This general mechanism arises in diverse branches 
of physics, such as gelation [Q, island formation in epi- 
taxial surface growth and stellar evolution By 
a long-term research effort, considerable understanding 
of this irreversible aggregation process has been achieved 
In this Letter, we focus on a different growth mech- 
anism that appears to provide a natural description for 
the evolution of city populations. This is preferential 
evaporation from smaller aggregates and preferential con- 
densation onto larger aggregates 

There are many examples where this evapora- 
tion/condensation mechanism occurs in physics and in 
the social sciences. The classic physics example is the 
late-stage coarsening of a binary mixture in an off-critical 
quench below but near the coexistence curve [^,0 . Here 
the system separates into droplets of the minority phase 
that are embedded in a matrix of the majority phase. 
Subsequent growth proceeds through preferential evap- 
oration from smaller droplets, due to the effect of sur- 
face tension, and subsequent condensation onto the larger 
droplets 

In the social sciences, it has been argued [|j that the 
growth of cities may be due to migration from small to 
large cities, as opposed to a view that emphasizes dif- 
ferential population growth ^. In a spirit closer to our 
work, the wealth distribution of individuals was described 
by a kinetic asset exchange model with preferential trans- 
fer from poor to rich individuals in each transaction [ p^ . 
For generic situations, the solution to the rate equation 
showed that this mechanism gives very different kinetic 
behavior from conventional aggregation. 

Motivated by these fragmentary results, we investigate 
a general class of migration-driven growth phenomena 
and show that, at large times and for large aggregate 
sizes, a comprehensive scaling theory can be developed 
with a minimum of assumptions. This theory gives both 



the growth rate of the typical aggregate size, as well as 
the asymptotics of the aggregate size distribution. Basic 
features of our theory agree with data on the popula- 
tion dynamics of U.S. cities. An important feature of 
our theory is that one may infer the general form of the 
migration rates from observations of the aggregate size 
distribution. For phenomena such as the city popula- 
tion distribution or the wealth distribution, we may thus 
hope to predict basic aspects of the dynamics in systems 
for which we have little a priori knowledge of underlying 
microscopic driving mechanisms. 

The model that we study is defined as follows. Let 
aggregates Aj be characterized only by their mass j, or 
equivalently, by the number of individuals that comprise 
them. These aggregates evolve according to the following 
irreversible reaction 



Ak + A, 
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That is, a monomer (or equivalently, one person) leaves 
a smaller aggregate of size k and joins a larger one of size 
I with rate K{k; I). This generalizes the asset exchange 
model of Ref. flQ], where a restricted class of reaction 
rates K{k; I) were considered. More generally, migra- 
tion could also go from a larger to a smaller aggregate. 
The symmetric limit, where the migration direction does 
not depend on the relative sizes of the two aggregates, 
leads to a diffusive-like kinetic universality class. We de- 
fer the investigation of this general system to a future 
work. Instead, we focus on the situation where there is 
preferential migration from small to large aggregates. In 
fact, any migration bias leads to scaling behavior for the 
aggregate size distribution identical to that of complete 
bias, as embodied by Eq. (Q). 

We now make the assumption of spatial homogeneity, 
so that the system is fully characterized by the concen- 
trations Cj{t) of aggregates Aj of size j at time t. We also 
assume that the law of mass-action applies so that the 
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time dependence of the aggregate concentrations may be 
described by the foUowing rate equations: 



oo 

)W = :j K{k;l)ck{t)ci{t) + <5(,j-i - <5fcj - (5; j] 

(2) 



The various delta-function terms enforce the constraint 
that the size of each aggregate changes by ±1 in a sin- 
gle reaction. The initial condition may be taken to be 



,(0) 



but any initial condition may equally well 



be considered, provided it is rapidly decaying in j. 

From these equations, we can immediately draw sev- 
eral important conclusions. First, there are no equilib- 
rium solutions. Rather, the size of aggregates grows con- 
tinuously and each Cj {t) eventually goes to zero as t cxd. 
Second, the total mass contained in the aggregates is (for- 
mally, at least) conserved. That is. 



(3) 



if the necessary interchanges between the infinite sums 
in this equation can be justified. Here we shall confine 
ourselves to this mass-conserving case. For definiteness 
we normalize the total mass to unity. 

We now make the conventional scaling ansatz for the 
large-time behavior of Cj(t) ^. We assume that there 
exists a well-defined typical aggregate size s(t) at time t 
such that 



(4) 



Here the exponent —2 follows directly from the condi- 
tion that the total mass is conserved, as discussed, for 
example, in [Q. We further assume that the reaction 
rates K{k; I) are homogeneous of degree A, or at least, 
that they are asymptotically so in the limit of large sizes. 
That is 



K{ak; al) = a^K{k; l)[l + o(l)]. 



(5) 



In the context of city population growth, the homo- 
geneity exponent A can be given the following interpre- 
tation. When the populations of two cities are scaled by 
some factor, there are both more susceptible migrants in 
the smaller city and potentially more reasons to move to 
the larger city. It is then natural that the overall migra- 
tion rate varies as a power law in this scale factor. Ex- 
ceptions to this behavior typically involve the existence 
of a cutoff size that separates two qualitatively different 
kinds of behaviors. 

Substituting the scahng ansatz (0) into Eq. (||), we 
find that s(t) satisfies s{t) = s{t)^ , with asymptotic 
solution, for A < 2, 



Defining z as the growth exponent of s(t), we thus have 

This growth exponent can also be obtained by adapta- 
tion of a back-of-the-envelope estimate for the typical 
size in irreversible aggregation. In aggregation, the re- 
action of aggregates of typical size s leads to a growth 
As of the order of s in a time At of the order of 
1/ (concentration x s^^). Here the concentration scales 
as 1/s and is the inverse reaction rate between 
typical-size aggregates. This leads to s ^ s'^, from which 
s ^ t'^/'-'^'^K For migration-driven growth. As is now of 
the order of 1 in the time At. This gives s ~ s^~^, thus 
reproducing the growth exponent of Eq. (0) . 

On the other hand, if A > 2, a power-law decay of the 
Cj(t) in J sets in at finite time. This feature invalidates 
the mass conservation statement and hence the scaling 
form of Eq. (||). The limiting case A = 2 can be treated 
within our scaling formulation, but must be handled with 
particular care, as we discuss below. This pattern of be- 
havior for the time dependence of the typical size paral- 
lels that of conventional aggregation, except that the size 
exponent in aggregation is z = 1/(1 — A) and a finite-time 
gelation transition occurs for A > 1 Note also that 

when the migration rate is symmetric, a scaling analy- 
sis similar in spirit to that just presented shows that the 
mean aggregate size grows as ^1/(3-^) for A < 3. Thus 
even migration without population bias leads to growing 
aggregates, albeit at a slower rate than if a bias towards 
larger aggregates exists. 

Also from the scaling ansatz, we find, after some non- 
trivial algebra, that the scaled aggregate size distribution 
$(x) obeys 



d<^>{x) 

dx 
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[Kix;y)-K{y;x)]^y). 



(8a) 



(8b) 



From these equations the basic qualitative behavior of 
^(x) can be deduced. Note that only the antisymmetric 
part of K{x; y) contributes to the scaling limit. Without 
loss of generality, we can now assume that K{k;l) — for 
k > I. From Eq. (|8^) it follows that $(a;) can be discon- 
tinuous whenever ^!{x) + x becomes zero. In particular, 
at such a point Xc, $(a;) can be consistently set to zero 
for all X > Xc- As long as $(a;) is different from zero, 
however, (pa) can be integrated to yield 
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(9) 



sit) cx [(2- A)t]'/^'"^^ 



(6) Here A and xq are arbitrary constants chosen so that 
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An important feature of <i>(j:) is its behavior for small 
values of a;. To quantify this, we define the exponent r 
through 



$(a;) cx x^-^ [1 + 0(1)] 



(11) 



With this definition, one has Cj{t) ^ j for 1 <C :/ <C 
s{t), as well as 



— (2 — t)z — 



(l«j«s(t)). (12) 



This defines the exponent w. 

To proceed further, we introduce another fundamental 
exponent that completes the scaling characterization of 
the reaction rates K{k;l), namely. 



K{1;1) Kl^-i" {I ^00). 



(13) 



This is entirely analogous to the corresponding defini- 
tion in conventional aggregation where the form of the 
cluster size distribution depends on the relative rates 
of small-small, large-large, and large-small reactions [||. 
With these definitions, we find, after detailed analysis of 
Eq. (^, four different classes of behavior: 

• Type 1: A>l,/i>l. In this case t — X and 
hence w = 1. 

• Marginal: X > 1, ^ — 1. Here it is not possible 
to make simple statements about the value of r. 
Rather, r depends on the complete shape of ^{x) 
and therefore on the very specific form of the reac- 
tion rates. 

• Type 2a: A < 1, /i < (1 + A)/2. In this case r = ^ 
and u; = (2 - /i)/(2 - A). 

• Type 2b: A < 1, (1 + A)/2 < ^ < 1. In this case 
T = (1 + A)/2 and w ^ {3 - A)/(4 - 2A). 

For the complementary large-x behavior of $(a:), we 
now show that in almost all cases $(x) vanishes beyond 
a certain critical value Xc of its argument. Indeed, sup- 
pose the contrary. It then follows from Eq. (|^) that 
\E'(x)^— cxDasa;— >oo. This can happen in three ways: 
either ^'(a;) varies faster than linear, slower than linear, 
or linearly in x. In the first case, Eq. (|^) would indicate 
that $(a;) < for large x, which is impossible. In the 
second case, $(x) would go to a constant as a; ^ 00, in 
contradiction to Eq. (p^. Thus the only viable possibil- 
ity is the third case, which occurs if X — fi = 1. 

A more thorough investigation is required to determine 
whether it is possible to find a consistent large- a; behav- 
ior for $(a;) in this last case. If so, then ^{x) would 
have a power-law decay such that the integral in Eq. ( |lO| ) 
still converges. However, in all other cases, Eq. (^ must 
cease to be valid at some point and the function ^(x) 



must vanish identically afterwards. This can happen in 
two different ways: Either x + ^(x) has a simple zero at 
some point Xc and the function ^{x) jumps from its value 
at Xc to zero, which is possible according to Eq. (|8a|), or 
else the function ^(x) goes smoothly to zero at Xc as a 
consequence of Eq. (^ and the double zero of a; + ^'(a;) 
at Xc- These results closely correspond to those of the 
Lifshitz-Slyozov- Wagner (LSW) theory of coarsening 
as well as to models of asset exchange . 

A special case that can be solved exactly is the case 
^ = X. For this situation, we find 



2-A 



{X < (2 

otherwise. 
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l/(2-^)^ 
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If A > 1, this case belongs to a system of Type 1 listed 
above, whereas if A < 1, this case belongs to Type 2a. For 
either alternative, the correct exponent r is predicted. 
Note further the discontinuity in ^{x) that indeed oc- 
curs exactly at the point where x + '^{x) vanishes. Many 
other cases can be handled similarly and will be presented 
in a forthcoming publication |ll| ]. 

A situation that requires a more refined analysis is 
A = 2 and /i > 1. For these parameter values, it follows 
that T = 2, which is incompatible with the normalization 
condition Eq. ([To|). To obtain valid results, we need to 
modify the scaling ansatz as follows 



Ins(t) 



It follows that 



(t) = exp y%t + B] 



(15) 



(16) 



where B is some constant. The function ^(x) then has 
the normalization $(0) — 1 and satisfies a modified ver- 
sion of Eq. (pa|). 

Let us now discuss how our scaling theory applies to 
LSW coarsening For this system, the migration rate 
K{i;j) is given by the product of the rate at which a 
particle evaporates from an aggregate of mass i and the 
probability that it reaches an aggregate of size j. In 
the evaporation step, the diffusive current J is Ac/R(i), 
where Ac is the difference between the monomer concen- 
tration near the interface and in the bulk. This difference 
is proportional to R{i)~^, since it is due to surface ten- 
sion. The current J is therefore of the order of R{i)~'^ 
and thus the rate at which particles leave an aggregate 
of size i is proportional to JR{iY~^, that is, to R{i)'^~^ . 
Further, the probability of reaching an aggregate of size j 
in three dimensions is simply proportional to its volume 
R{jY- We therefore find for the overall migration rate 



X(z;j)«i?(z)''-3i?(j)'^«z(^-3)/^j-. 



(17) 



From the definitions of A and fi in Eqs. (|) and (|l|), 
it follows that the system is of Type 2a, from which one 
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obtains z = d/3 and r = 1 — 3/d. These indeed corre- 
spond to the LSW predictions, in which the characteristic 
cluster radius increases as t^^^ Q and the number n{R) 
of clusters of radius R ^ Rc{t) varies as R^ Our 
theory also correctly predicts that the scaling function 
vanishes beyond a certain value of the scaling variable. 
On the other hand, the migration rate of Eq. is not 
precise enough to ensure that our scaling theory repro- 
duces the same functional form for $(a;) as that of the 
detailed LSW theory. 

To apply our theory meaningfully to the evolution of 
city populations, it is necessary to incorporate the effects 
of demographic population growth. Over intermediate 
time scales (of the order of decades) , demographic growth 
typically gives a population that increases exponentially 
with time. Such a behavior can be modeled by allow- 
ing the process — > A^+i to occur at rate kj. When 
demographic growth and migration occur together, the 
scaling ansatz for the underlying rate equations needs to 
be modified accordingly. We have found that the appro- 
priate scaling ansatz for this more realistic situation is 



(18) 



With this hypothesis, the functional form of <i>(x) turns 
out to remain the same as the case of no demographic 
growth, but now the typical city population grows as 



„7t/(2-A) 



(19) 



Hence we arrive at the central conclusion that the typ- 
ical city size grows much faster than the population of 
the country as a whole as A approaches 2. 

The city population distribution in many countries is 
consistent with a power-law form in which the exponent 
T is close to 2 Our scaling theory then requires 

that the homogeneity exponent A is also close to 2. Thus 
from Eq. (p^), the typical city population should increase 
much faster than the overall population. This is con- 
firmed qualitatively by data for the populations of various 
U.S. cities during their early histories ||l^. The popula- 
tion of essentially every major U.S. city grows much faster 
than the U.S. as a whole over considerable time range. 
However, as cities reach maturity, their growth may slow 
or their population may even decline for reasons unre- 
lated to preferential migration to still larger cities. 

In summary, we have introduced a simple kinetic de- 
scription for migration-driven growth and developed a 
scaling theory that determines the large-time behavior 
for the aggregate size distribution. Asymptotic results 
depend only on rudimentary properties of the reaction 
rates, most notably the homogeneity index A. The typi- 
cal aggregate size grows as t^^^^^^\ while several distinct 



behaviors emerge for aggregate size distribution. Our re- 
sults represent the counterpart of the scaling theory of 
irreversible aggregation to migration-driven growth. Fi- 
nally, we have suggested a connection between migration- 
driven growth to the distribution of city populations and 
found a qualitative correspondence between model pre- 
dictions and recent data on U.S. cities. 
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DMR9978902, and DGAPA INI 12200 for financial sup- 
port of this work. 
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